Optimized Dynamical Decoupling for Power Law Noise Spectra 
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We analyze the suppression of decoherence by means of dynamical decoupling in the pure- 
dephasing spin-boson model for baths with power law spectra. The sequence of ideal tt pulses 
is optimized according to the power of the bath. We expand the decoherence function and sep- 
arate the cancelling divergences from the relevant terms. The proposed sequence is chosen to be 
the one minimizing the decoherence function. By construction, it provides the best performance. 
We analytically derive the conditions that must be satisfied. The resulting equations are solved 
numerically. The solutions are very close to the Carr-Purcell-Meiboom-Gill (CPMG) sequence for 
a soft cutoff of the bath while they approach the Uhrig dynamical- decoupling (UDD) sequence as 
the cutoff becomes harder. 

PACS numbers: 82.56Jn, 82.56.Dj, 03.67.Lx, 03.65.Yz 



I. INTRODUCTION 

The dynamics of a spin S = 1/2, or a quantum bit 
(qubit), coupled to an environment is one of the longest 
studied models of quantum decoherence. It finds im- 
portant applications both in nuclear magnetic resonance 
(NMR) and in quantum information processing (QIP). 
The suppression of the decoherence is one of the goals 
one usually strives for. In order to achieve it, sequences 
of control pulses are used 

There are two ways to address dynamical decoupling 
(DD) by control pulses. Either the control is modulated 
continuously, see for instance Ref. @, or the control con- 
sists of short pulses which can be seen as (approximately) 
instantaneous. In this article, we will focus exclusively 
on the latter approach which is common in NMR and 
wide-spread in QIP It relies on control pulses to 

invert the dynamics of the spin. The spin experiences a 
rotation of an angle tt at each pulse. The choice of the ap- 
propriate sequence is essential for an enhancement of the 
suppression of decoherence. A large variety of sequences 
has been suggested. The majority is characterized by 
periodic pulses [1, 0, H, 0,11, Q- The most famous 
example is the Carr-Purcell-Meiboom-Gill (CPMG) se- 
quence d, [§| where cycles of two n pulses are iterated. 

Other sequences proposed are non-equidistant such as 



the concatenated dynamical decouphng (CDD) 17, jOl- 111 
or the Uhrig dynamical decoupling (UDD) [ij, [H, fl4| . 
The CDD consists of concatenations of pulse sequences. 
It can suppress both transverse relaxation and longitu- 
dinal relaxation at the price of a relatively large number 
of pulses. If £ is the largest order in an expansion in 
the total duration t, in which no decoherence occurs, the 
required number of pulses grows exponentially with i. 
The UDD eliminates only pure dephasing, but in turn it 



requires only a linearly growing number of pulses. The 
concatenation of the UDD sequence (CUDD) allows for 
the suppression of transverse and longitudinal relaxation 
again at the price of an exponential growth of the num- 
ber of pulses, but requiring only the square root of the 
number of pulses necessary for CDD |14l |. 

We point out that all the above sequences are idealized 
in the sense that they are based on ideal, instantaneous 
pulses, i.e., S peaks, though the effect of finite pulse du- 
rations is being discussed. Furthermore, sequences of re- 
alistic pulses have been proposed [H, [H, [o, (HI . The 
finite duration of the pulse is a source of additional er- 
rors which can be reduced by designing the shape of the 
pulses appropriately, see for instance Ref. [l^ and refer- 
ences therein. In the present paper, however, we will 
concentrate only on sequences of instantaneous tt pulses. 

The UDD sequence was discovered first for a spin- 
boson model where it was observed that no details 
of the model entered. On the basis of numerical evidence 
and finite order recursion it was co nject ured that UDD 
is applicable to any dephasing model [l3l . [20l | . This claim 
was finally proven for arbitrary number of pulses in the 
total duration t of the sequence 2l|. For various clas- 
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sical noise spectra the experimental verification of the 
theoretical results was achieved by optical control of the 
transition in Be ions [13, HI, [13| • It was also shown that 
the UDD sequence outperforms the CPMG sequence and 
equidistant sequences in general for pure dephasing baths 
with hard cutoff while it performs worse for soft cutoffs 
[H HO, mill Hi]. This shows that the knowledge of the 
cutoff is an essential piece of information for an optimum 
suppression of decoherence. 

The performance of some pulse sequences for classical 
noise spectra has already been considered by Cywihski et 
al. for superconducting qubits subject to gaussian and 
random telegraph noise ^25]. The authors compare the 
efficiency of different sequences in suppressing pure de- 
phasing. They find that the UDD sequence is optimum in 
suppressing the decoherence if the gaussian noise displays 
a hard ultraviolet (UV) cutoff. In situations, however, 
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where one has to work in the regime of smah frequencies 
(long times) such that the cutoff cannot be reached, the 
CPMG sequence is the one yielding the best results. 

Another proposal for an optimized sequence is put for- 
ward by Biercuk et al. [2^ l23l | . The optimization in the 
UDD is extended to a locally optimized dynamical de- 
coupling (LODD) sequence which is tailored to a given 
experimental noise environment. The experimental im- 
plementation of LODD for classical noise shows that it 
performs better than UDD and CPMG. The LODD, how- 
ever, is limited by the degree to which the spectral func- 
tion of the noise (or the bath) is known. This caveat 
is dealt with by an optimized noise-filtration dynamic 
decoupling (OFDD) [2J| where the noise spectrum is ap- 
proximated by a constant up to a high-energy (UV) cut- 
off for which sequences of pulses are deduced numerically. 
For an ohmic bath the OFDD implies about the same fac- 
tor of w 1.5 of improvement over the UDD sequence than 
the more cumbersome LODD. But it does not provide a 
significant improvement for an ambient noise scaling such 
as oc 

In this paper we generalize the OFDD from a constant 
spectrum to an arbitrary power law without UV cutoff. 
This problem can be analysed to a large extent analyt- 
ically. The resulting optimized dynamic decoupling se- 
quences for power law spectra (PLODD) are universal 
in the sense that the relative instants {6j = tj/t} of the 
pulses depend only on the power law exponent a. The 
gist of our finding is that the PLODD sequences resem- 
ble CPMG sequences for slowly decreasing noise spectra 
while they approach UDD sequences for fast decreasing 
noise spectra. This agrees with the qualitative expecta- 
tions based on other investigations l22|, ^2^, ^2^, i2S\ . 
Our results provide an important guideline in the choice 
of sequences to be applied in experiments on the suppres- 
sion of decoherence. 

From an experimental point of view, the noise used to 
test the dynamical decoupling sequences - classical noise, 
ohmic spectrum, as well as the l/w^ noise spectrum - 
is in general governed by a power law spectrum. Ana- 
lytically, power law spectra approximate any spectrum 
either for very small or for very large frequencies which 
in turn correspond to long or short durations, respec- 
tively. Moreover, power law spectra are the perfect tool 
to analyze the influence of baths characterized by differ- 
ent cutoffs on the optimized pulse sequence. By varying 
the exponent of the power law of the spectrum we can 
simulate both baths with soft and hard cutoffs and we 
can interpolate smoothly between them. 

In order to render an analytical investigation possible 
we consider the spin-boson model with pure dephasing. 
We start from the decoherence function xit) (defined in 
Eq. (HI)) which measures the size of the decoherence. The 
advantage is that the decoherence function at the same 
time embodies both the characteristics of the bath and 
of the sequence of tt pulses. Hence one has to minimize 

xit). 

The article is organized as follows. In Sect. |TT] the 



model is introduced and the decoherence function for a 
sequence of tt pulses is defined and discussed. In the fol- 
lowing Sect, mil we study and solve the problem of the 
diverging terms in the integration yielding the decoher- 
ence function. Then we derive the final equations in Sect. 
IIVI and solve them numerically in Sect. |Vl where also ex- 
amples of the PLODD are shown. At last the conclusions 
are drawn in Sect. |Vll 

II. CONVERGENCE OF THE DECOHERENCE 
FUNCTION 

We consider the spin-boson model with pure dephasing 

H^Y. + E ^' +b^)+E (1) 

i i 

describing a single qubit as a spin 5=1/2 coupled lin- 
early to a bosonic bath. The spin is represented by the 
Pauli matrix cr^, while the b[^^ are the annihilation (cre- 
ation) operators of the bath. The constant E sets the 
energy offset. The properties of the bath are defined by 
the set of parameters {A^, w^}. This information is con- 
veniently encoded in the spectral density [H, [l^] 

J(c.)-^A2^(c.-c.O. (2) 

i 

We recall that the quantum mechanical time evolution 
p" of a sequence with n tt pulses about the x axis of the 
spin reads 

P" = /t-5„^/<5„-5„_i^- • ■ X fss-S2X fs2-5iX fs^, (3) 

where X stands for the spin operator of the rotation due 
to the pulse [13, [21I . Such a sequence suppresses the 
relaxation along z jl4l.[2l|. If i is the total duration of the 
sequence, the instant tj, at which the pulse j is applied, 
is given by tj = tSj. By definition Sq := and ^n+i ■= 1 
although there is no pulse neither at the very beginning 
nor at the very end. The notation fsi-St^x stands for the 
free evolution of the system in the interval t{5i — 5i-i) 
between two successive pulses. 

In Refs.[l^.[T3lit is shown that the free induction decay 
is proportional to e~^-^^*^ where the decoherence function 
is defined by 

X(0:= r^lynMPdc.- (4) 

Here the noise spectrum S{lij) is related to the spectral 
density J(a;) in by 

5(w) := ij(u;)coth(/3w/2), (5) 

where /3 is the inverse temperature. The filter function 
?/„(z) (z := Lot) for n pulses is given by 

n+l 

y„(z):=^2«^■(-l)^"e-^^ (6) 
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with 



if j 0, n + 1 

1 if j e {1,2,.. .,n- l,n} 



(7) 



Obviously, it encodes the properties of the sequence. 

Equation ([4]) is the starting point for the evaluation 
of the optimized sequences. The aim is to keep e"^-^^*-* 
close to the unity as long as possible. In Ref. [13 the 
condition was enforced that the first n derivatives of the 
filter function should vanish a.t ujt — for a sequence 
with n pulses. Thus the function y„ would increase very 
slowly close to zero. The condition on the derivatives 
implies the following set of non-linear equations 



n+l 



(8) 



for p e {1, 2, n}. For p = 0, 
which is equivalent to j/„(0) = 
Eqs. © reads [l3| 



Eq. (O is also zero 
0. The solution of the 



rUDD 



[j7r/2(n + 1)] . 



(9) 



The UDD sequence suppresses the decoherence best for 
baths with a hard cutoff rather than for baths with a very 
soft cutoff. This was tested by means of numerical simu- 
lations [lOl , analytical analyses [l^, [25| , and experiments 
with classical noise [22, . 
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Figure 1; (color online) Schematic visualization of the filter 
function (black dashed line) and of the spectral function S{iij) 
for a soft and for a hard cutoff, ujd stands for the ultraviolet 
cutoff. The plot illustrates that the overlap between the filter 
function and the spectral density remains significant for a soft 
cutoff, even if |j/„([i;t)p shifts to the right. This is in contrast 
to the situation for a hard cutoff where the overlap becomes 
extremely small for n cx> or t 0. 

The efficiency of UDD depends on the applicability 
of an expansion in powers of t of x(t) [H, [23|- This 
expansion implies the expansion of yn{z) in powers of 
z = ujt. This is always possible since yn{z) as defined 



in ^ is analytical. But the existence of the integrals of 
the resulting series in powers of wt, as required by ([3]), 
depends on the UV cutoff of S{ui) 

In the derivation of Eq. ^ only the existence of the 
derivatives of ?;„ {ujt) and not the existence of the integral 
over the frequency uj is required. The decoherence func- 
tion consists of the product of the function S{ijj) times 
the square modulus of the filter function. The function 
x{t) is minimum if the overlap between S{uj) and 
is minimum. The significance of the UV cutoff is illus- 
trated in Fig. [1] 

From now on we focus on power law spectra 



S{uj) So 



,a+i ■ 



Then the decoherence function reads 



1 



— |y„(tji)| dw 



(10) 



(11) 



It is required that a is strictly positive to ensure the 
convergence for uj —^ oo. This is true because the filter 
function is bounded from above \yn{ujt)\ ^ 2(n-|- 1). The 
prefactor 5*0 incorporates all the constants of the spectral 
density. 

For Lu ^ (infrared (IR) limit), x(^) converges if 
yn{ujt) oc (ujt)"^ for m large enough. This in turn de- 
pends on the choice of the sequence {6j}. For arbitrary 
sequences {5j} we have yn{0) = so that the IR conver- 
gence is guaranteed for a < 2. For larger a, we require 
that the first m derivatives of the filter function vanish, 
implying 



|y„(c.t)pa(c.t)2(™+i). 



(12) 



which is similar, but not identical, to the requirement for 
the UDD in Eq. ([8|). The convergence of x(t) is guaran- 
teed for 



a < 2m 



(13) 



In Ref. it was argued that UDD applies indepen- 
dently of the temperature. This indicates that UDD 
can be equally used to suppress classical gaussian noise 
[H, • For high temperature (3^0 the thermal fluctu- 
ations dominate over the quantum fluctuations such that 
S{uj) (X l/u! for J(0) 7^ 0. This is the famous 1/f noise, 
which corresponds to a = 2 in our notation. The case 
S{uj) oc l/oj'^ is experimentally relevant for ions in a Pen- 
ning trap [23, HI] . This case corresponds to a = 5 in the 
above notation. 

The basic 2-pulse cycle of the CPMG sequence coin- 
cides with the UDD sequence for n = 2. It makes the 
first two derivatives of the filter function vanish [l2| . Ac- 
cording to Eq. [T3] its applicability is restricted to baths 
characterized by a < 6. 

In order to study general power laws we proceed as 
follows. We substitute z = ujtin the decoherence function 
x{t) in ^ obtaining 



xit) = Sorin 



(14) 
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with 



In 



dz. 



(15) 



This simple substitution reveals that the optimum {Sj} 
are independent of the total duration t of the sequence. 
All the time-dependence of x{t) is a simple power of t 
as in (|14p . Its exponent is determined by the power law 
spectrum of the bath. The precise sequence {Si} deter- 
mines the factor /„. 

The condition for the first m derivatives to vanish is 
given by the set of the first m non-linear equations in 
([5]), i.e., for p e {1,2,..., m}. These conditions are the 



same as those leading to the UDD sequence [12| except 
that they do not need to be fulfilled up to p = n but only 
up to p = m. For a sequence of n > m pulses we still 
have n — m degrees of freedom left. This freedom is used 
to minimize /„ and hence the decoherence function x(t). 
In this minimization the m Eqs. ([8]) act as additional 
constraints. Hence, we have to study the variation 



_d_ 



In 



dVn{z) 



4=1 



dz 



(16) 



The regularization by a finite IR cutoff x is required 
because the limit a: — > does not exist for the individual 
terms lij {x) . Each term [x) can be reduced to an 
analytical expression by the substitution z — > —z//S.ij 



h,{x) = (-A,,r 



dz 



= (-Ay-)"r(-a,-A„a;) 



(20a) 
(20b) 



where r(— a, — Ayx) is the incomplete Gamma function 

For later use we state that the vanishing of the first m 
derivatives of the filter functions yn{z) at z — implies 
(fn)) and thus 



n+l 

E 

ij=0 







21^+1, {-iy+3 (Ay f (21a) 
(21b) 



for ^ p ^ 2m + 1. Eq. (|21b[) asserts that the weighted 
sum of powers of A^ vanishes. We will utilize this can- 
cellation to find the relevant contributions to /„ analyti- 
cally. 



where m Lagrange multipliers Xi appear due to the m 
constraints. 



A. Expansion of the incomplete Gamma function 



III. DIVERGING TERMS 

The integral /„ (fT5|) converges if the condition a < 
2(m-|- 1) is fulfilled. But the integration in ^ cannot be 
carried out analytically. To make analytical progress we 
split the square modulus of the filter function into a sum 
of exponential terms according to 



n+l 

\yn{z)f= E 2' 



(-1) 



where we use the notation Aij :— i{Si — Sj) — 0+ for 
{i, j} = 0..n -I- 1; the summand — 0"*" stands for an in- 
finitesimal negative real part which is required later on 
for convergence for z +oo. The integral /„ is given by 
the limit of the sum 



/„ = hm In{x) 

n+l 



(18a) 



where the integrals 



Iij{x) := / —T-rdz. 



(19) 



UV convergence is ensured by the infinitesimal negative 
real part of Ay , see definition below Eq. P7|) . 



a. Non-Integer Exponents Starting from the defi- 
nition of the incomplete Gamma function r(a, x) with 
—a = a ^ No we integrate by parts £ times and write 



r(a,x) = Ti + 7^. 



with 



e-i 



p=0 



a+p (a - 1)! 
{a+p)l 



(17) and 



To 



(g^ 1)! r{a + £,x) 

(a + i-iy. ■ 



(22) 



(23) 



(24) 



For a ^ the limit lim^r^^o r(a, x) is not defined. We 
choose £ such that 



a + £>0> a + £-l. 



(25) 



For non-integer a this is possible. The first inequality 
ensures that limx^o T{a+£, x) = T{a+£). The recurrence 
relation T{z -|- 1) = zT{z) implies 



lim T2 — r(a). 



(26) 



Now we concentrate on the term 7i. Since the exponen- 
tial function can be expanded in powers of x, Ti has 
a well-defined expansion 



Ti{x)=x''"YapxP. 

p=0 



(27) 
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The coefRcients {ap} depend on the coefRcients of Ti 
and on those of the expansion of the exponential. Their 
exphcit form does not matter here. The powers diverge 
in the hmit x for p £ {0, ...,£— 1} . For p ^ £ they 
vanish for a; — > 0. 

Finally we write the integral In{x) in terms of the in- 
tegrals Iij{x) defined in Eq. and evaluated in Eq. 
(l2Ub)l 



i-i 

In{x) = a^x 

p=0 
n+l 



n+l 

,p-a 2'?'+«^(-l)*+^'(-Ay)P (28a) 



+ \ 2'-+«^(-l)^+-'-(-A.,rr(-a) 



(28b) 



+ ©(x^-"). (28c) 

The last term ©(x^^") comprises all the contributions 
which vanish for a; — > 0, recall (P5|) for a = —a. The 
first term (|28ap contains the diverging terms. But they 
cancel one another completely because the inner sum in 
P8ap vanishes due to (j21bp . To sec this one must use 
(fni and ^ to arrive at p < ^ - 1 < 2m + 1. This 
rather formal argument simply reflects the fact that (fT51) 
guarantees the IR convergence of the integration in 
Hence all IR divergent terms appearing in intermediate 
calculations have to cancel finally. 

From the above the only remaining and thus relevant 
contribution to /„ = lima;_>o+ In{x) is (j28bp 

n+l 

/„ = r(-a) 2''+'^(-l)*+^(-A„r. (29) 

This is the result for non-integer a. Note that T{—a) is 
only a global prefactor which does not depend on A^- . 

h. Integer Exponents Next we consider the case of 
a = No. We use Eqs. (|23l24p with I = a = -a 



T{-a,x)= - e-"|]xP-°(-l)^+P ^" \ 
(-1)" 



a 



-r(o,x) 



(30) 



Thus Eq. ([28]) now reads 

n+l 



a-1 

In{x) = ^ apX' 
p=0 



^ 29-+«^(-l)*+^(-Ay)f (31a) 



where 7 is the Euler-Mascheroni constant. Hence in the 
limit a; — s- the only non-vanishing contribution to /„ = 
limj;^o-^n(a;) is 



n+l 



^" = 7;r E 2^'+''^(-l)^+-'Ar^.ln(-A.,). (33) 



ij=0 



B. Example for a = 2 



For a = 2 we explicitly write the expansion of the 
incomplete Gamma function and show that the diverging 
terms cancel. For a; — > we have 

A2. 

A2^T(-2,A,,a;) = ^(3 - 27 - 2 ln(-A,,a;)) 

+ ^ + ^ + (34a) 

A2. 

__^ln(-Ay) + 0(a;). (34b) 



The function li''\x) is the sum of the I^j'ix) weighted 
according to the right hand side of (jl8bp . These contri- 

( — 2) 

butions vanish. For /„ we can write 



(fe)/ 



i,j=0 
n+l 



-J (fc-Tj. n+l 



2a;2 



i=0 



3=0 



= 0. 



(35c) 

The vanishing of P5b|) is guaranteed by Eq. ^ for p — 
or, equivalently, by the property y„(0) = 0. Because of 
the sum over all {ij} we have 



n+l 



^n""(^) = - E 2^'+«^ (-l)^+^(<5. - 6,) (36a) 

n+l n+l 

^2^^(-l)^-^2''(-l)M. 



I 

X 



3=0 

n+l 



1=0 
n+l 



_^2^'(-l)'^2'^(-l)^<5, (36b) 
0. (36c) 



a! 



2?»+9. (-1)*+J(-Ay )"r(0, -Aya:).(31b) Similarly one obtains for 4/ 



(0) 



i,j=0 



As before the inner sum in (|31ap cancels because of (|21bp 
and because (fT5)l implies p < a — 1 < 2r7i. Note that (fT^]) 
and (|21bp additionally imply that the weighted sum of 
the powers Ag vanishes. This will help to simplify (|31bp 
further in the limit a; — s- 0. We expand the incomplete 
Gamma function 

r(0, -A„-x) = -7 - H^) - ln(-A,y) + 0{x), (32) 



27 + 21n(a;) -3 



n+l 



^2«-+«^(-l)'+-'"(<5,-J,)' 

= 0. (37) 

The last equation vanishes because of Eq. ([8]) at p = or 
at p — 1. This is because (S^ - 5jY = Sf + - 2S.iSj. Eq. 
(ISl) has to hold for p = and p = 1 because these values 
fulfill p < m since m > 1 is required by p3p for a = 2. 
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IV. RELEVANT TERMS AND FINAL 
EQUATIONS 



Eqs. (|29|33|) provide the analytical results for /„. One 
further simplification stems from the fact that /„ is real. 
Hence only the real parts of the summands in (|29l33p 
need to be included since the imaginary parts cancel. 
We analyze the case a £ No, where we distinguish even 
and odd exponents, and the case a ^ No separately. 

A. Integer Exponent 

c. a even. We use Ay :— iipij^O^ with ipij — Si—Sj 
and ln(re*^) = ln(r) +26* with r > and 16*1 < tt to obtain 

TT 



ReA^H-A.,) - (-l)"/VSRe(ln| 

^ (-1)"/V51n|(p.,|. 
The sum over all i ^ j as required by (|33p reads 



(38) 



reven _ 

n+1 



X 



^ 2'^^+'i^{-iy+^\5,-5,\-\r,\5,~5,\. (39) 



i,3=0 

d. a odd. Starting from 
ReA^ ln(-A,,) = 

{-l)^vf^Re[t (in 
= (-l)^|^.,r|. 
the final integral /„ becomes 



(40) 



odd 



2 



al 

n+l 



(41) 



i,j=0 




c^oPLODD 
X X CPMG 

UDD 
A ACDD 



10 12 14 



Figure 2: (Color online) Comparison of the optimized se- 
quences PLODD and UDD and the sequence CPMG for a = 4 
and for various numbers of pulses n. For clarity only the 
first half of the sequences is plotted but the symmetry about 
S — 0.5 (Sn+i-j — 1 — Sj) allows for the straightforward re- 
construction of the second half. The concatenated sequence 
CDD is shown only for n = 10 because it does not exist for 
n = 20 or n = 30. Note that we refer here to the CDD for 
pure dephasing fH^ . 
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B. a Non-Integer 

For positive non-integer a we consider 

Re (-Ay)" = |(^y|"Re (e-''(''/2)a sgn.^.,^ (42^) 
= cos((7r/2)a)|(^yr, (42b) 

which implies 

C = cos((7r/2)a)r(-a) 

n+l 

X J2 2''*+«^ (-l)''+^'|(5, - ,5^1". (43) 
where the superscript 'ni' of stands for 'non-integer'. 



Figure 3: (Color online) Log-log plot of the integral In as 
a function of the total number of pulses n for a = 2 and 
a — 4. Different sequences are compared, including the CDD 
sequence for pure dephasing [ll[. The latter takes only the 
values n = {0, 1, 2, 5, 10, 21, 42, . . .}. 



V. NUMERICAL RESULTS 



The power law optimized dynamical decoupling 
(PLODD) is a bath-optimized sequence depending on the 
exponent a only. For various values of a, we numerically 
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solve the system of non-linear equations 

ri+l 

= ^29^ (-!)■'■ 

[a/2] 



= 



d_ 



dVn{z) 



dz 



(44a) 
(44b) 



where p S {I,---, [§]}; here [x] stands for the largest 
integer not larger than x £ R and [a/2] results from 
(fO)) . The prefactor /„ is computed analytically both in 
the case of integer and of non-integer values of a; it is 
given in Eqs. (|39l41l43p . 

In the sequel, we restrict ourselves to symmetric se- 
quences for two reasons. First, the main results in dy- 
namical decoupling are derived for symmetric sequences. 
Second, we searched for asymmetric optimized sequences 
for small number of pulses, but those sequences found did 
not perform better than the symmetric ones. Hence we 
focus on symmetric sequences fulfilling 5n+i-j = 1 ^ <5j- 
In this case one has to deal with a system of n/2 + [a/ 2] 
equations, solved for n/2 variables [Si] and [a/2] La- 
grange multipliers Am. For example, for a ~ A the system 
consists of n/2 -I- 2 equations. 

Figure [2] shows the resulting PLODD sequences for 
a = 4 for various number of pulses n. The PLODD 
sequences are very close to the CPMG ones. We recall 



that 5\ 



CPMG 



(2j — l)/2n. No relevant dependence of 



the PLODD instants as functions of the number of pulses 
n can be observed. The concatenated sequence CDD for 
pure dephasing is also shown for comparison for n = 10. 
We recall its recursion p^^^ — p^^ Xp^^ for n even 
while Pn+i = holds for n odd; po stands for 

the free evolution without pulse. 

In Fig. [3] the evaluation of the prefactor /„ shows that 
PLODD performs shghtly better than CPMG for a = 4 
while for a = 2 the data for PLODD and for CPMG 
coincide. It is interesting to notice how the performance 
changes with a. For a = 2 we see that CDD and UDD 
provide almost the same results while for a = 4 UDD 
performs better than CDD, though still outperformed by 
CPMG and PLODD. In addition. Fig. [3] indicates that 
In decreases if a increases. We will come back to this 
point below. 

The log-log plot in Fig. [3] shows that /„ scales like 
a power law in n. The regression ln/„ ~ ailnn -I- ap 
yields oq = -(2.33 ±0.01) and ai = -(3.041 ± 0.003) for 
PLODD and ao = -(2.231 ± 0.009) and ai = -(3.062 ± 
0.003) for CPMG. We can compare these results with Eq. 
(25) in Ref. [2^ derived by Cywihski et ai. Cywihski's 
formula was derived for a de facto infinite UV cutoff for 
1.5 ^ a ^ 2.5 (or 0.5 ^ acyw ^ 1-5 in the notation of 
Ref. [25!) . It shows that x{t) scales like n"~^. 

As a further check we calculated /„ versus n for CPMG 
for 1/f noise, i.e., a = 2. We find ai = -(0.9899±0.0008) 
and flo = —(0.199 ± 0.002), which is close to the cor- 
responding value InCi ~ —0.163 reported in Ref. [2^. 



Cywihski et al. concluded that in the range 0.5 ^ 
acyw ^1.5 CPMG is to be preferred over UDD for the 
prolongation of qubit coherence. The UDD outperforms 
the other sequences in the range where a finite UV cutoff 
makes itself felt. Our systematic minimization confirms 
the results by Cywihski et ai. It extends them by putting 
them on a systematic basis leading to the optimum power 
law dynamic decouling and because a larger range of ex- 
ponents is treated. 

If we consider higher values of a, the PLODD sequence 
approaches the UDD one. This is illustrated in Fig. 
m The instants of the PLODD sequence appear to be 
bounded from below by the instants of the CPMG se- 
quence and from above by the those of the UDD one. 
For a > 6 CPMG does not satisfy the condition (fT5|) any- 
more which is required for the convergence of x{t) so that 
no comparison to PLODD or UDD is possible. For small 
values of a, [a/2] is also small. This implies that the con- 
ditions (|44ap are less important than the conditions (|44b[) 
resulting from the minimization of the prefactor /„. This 
finding confirms what was already expected from studies 
on UDD [1^, [23, [2^ , namely that the suppression of 
decoherence is more efficient for baths with a hard cutoff. 
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Figure 4: (Color online) Evaluation of PLODD sequences for 
various exponents for n = 10. For clarity only the 5i from 
i = 1 to 5 are plotted. The second half of the sequence is 
symmetric to the first one = 1 — 5j). The PLODD 

approaches UDD for increasing a. 

As a increases the PLODD sequences tend to coincide 
with the UDD sequences. We computed the maximum 
difference between the {5i] of the PLODD and of the 
UDD as a function of a. The results are depicted in Fig. 
O Clearly, the maximum difference decreases as a in- 
creases. This supports that the PLODD tends to recover 
the UDD for large values of a. This can also be seen in 
Fig.m 

The inset of Fig. [5] shows that In /„ decreases linearly 
with a for fixed values of n. Combined with the power 
law scaling shown in Fig. [3] we deduce that 

/„ (X (C„,„/n)"-i (45) 
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0.08 rr 




Figure 5; (Color online) The maximum difference between the 
{Si} for PLODD and UDD is shown for two sequences with 
n = 10 and n — 12 pulses. The PLODD sequences tend to 
approach the UDD sequences for harder cutoffs. In the inset 
the dependence of /„ on a is depicted. Note that the y axis 
of the inset is in logarithmic scale. 

holds for the PLODD sequences where the factor Cn,a 
depends only very weakly on n and a. 

VI. CONCLUSIONS 

We analyzed the supression of decoherence by means 
of sequences of instantaneous tt pulses. The work horse 
is the spin-boson model with pure dephasing which can 
be treated analytically. The sequences are optimized for 
power law spectra. The main difference to the already 
known optimized sequences, for instance UDD or OFDD, 
is that the convergence of the decoherence function x{t) is 
investigated for various powers (a) of the noise spectrum. 
Our approach extends previous results in several ways. 

The OFDD sequences proposed by Uys et al. are 
optimized for a constant spectrum (a = 1 in our nota- 
tion) with finite UV cutoff In our study we sent 
the UV cutoff to infinity and treated general power law 



spectra characterized by the exponent a > 0. Hence the 
proposed PLODD sequences are optimized for arbitrary, 
but fixed, exponent. They are universal in the sense that 
their relative switching instants Sj — tj/t depend only on 
a, but not on the total duration t of the sequence. 

In their investigation of classical noise with exponents 
a G {1.5,2.5}, Cywihski et al. observed that the well- 
known CPMG sequence works well in the regime where 
the UV cutoff is infinite for practical purposes. The UDD 
does not provide an improvement [1^. One of the authors 
generalized this investigation to the quantum mechanical 
spin-boson model and cutoffs with arbitrary power law 
behavior (Tsj . 

In the present work we extended these findings further 
by the systematic optimization of the sequence on the 
basis of analytical results for arbitrary power law spectra. 
The softer the UV behavior of the power law spectrum is, 
i.e., the smaller its exponent a is, the more the PLODD 
approaches the CPMG. Vice versa, the harder the UV 
behavior of the power law spectrum is, i.e., the larger 
its exponent a is, the more the PLODD approaches the 
UDD. 

Hence, the findings of previous investigations are cor- 
roborated. There is no completely different sequence 
which displays a significantly better performance for pure 
dephasing other than PLODD. The PLODD has the lim- 
iting cases CPMG (soft UV behavior) and UDD (hard 
UV behavior) . 

We also investigated how the decoherence function x(t) 
scales with a and n, the number of pulses. From Eqs. (|14p 
and (|45p we obtain 

X{t) cc n{Cn,^t/nr, (46) 

which generalizes the result in Ref. [2^ to arbitrary expo- 
nent and a quantum mechanical model. Here Cn,a is a 
factor which depends only weakly on n and a. 

We reckon that the results of this work represent a 
useful contribution to the technique of dynamic decou- 
pling. The optimization of the pulse sequence in relation 
to the specific baths is of vital importance in many ap- 
plications in high precision nuclear magnetic resonance 
and quantum information processing. 
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